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Exclusive Radiative Decays of B Mesons* 

Stefan W. Bosch'^ 

'^Max-Planck-Institut fiir Physik, Werner-Heisenberg-Institut, 
Fohringer Ring 6, D-80805 Munich, Germany 

We present within the Standard Model the exclusive radiative decays B — * K'/py and -Bs/d 77 in QCD 
factorization based on the heavy-quark limit rrib ^ Aqcd ■ For the decays with a vector meson in the final state 
we give results complete to next-to-leading order in QCD. 



1. Introduction 

The main goal of todays B physics is a pre- 
cision study of the flavour sector and the phe- 
nomenon of CP violation that comes along with 
it. Particularly useful for clean tests of the Stan- 
dard Model are the rare radiative b — > s(d)'y tran- 
sitions. The inclusive b ^ s'y mode was mea- 
sured to have a branching fraction oi B{B 
Xs-f) = (3.23 ± 0.42) • 10-*. Branching ratios 
of exclusive radiative channels are measured to 
be B{B° -> K*°-f) = (4.44 ± 0.35) • IQ-'^ and 
B{B+ K*+-f) = (3.82 ± 0.47) • IQ-'^ whereas 
for the B fTf and B — > 77 decays so far only 
upper limits exist 

Whereas the inclusive mode can be computed 
perturbatively using the heavy-quark expansion, 
for the exclusive channels bound state effects 
have to be taken into account. The basic mech- 
anisms for the exclusive radiative decays were 
already discussed by various groups [^. How- 
ever, they all had to use hadronic models which 
do not allow a clear separation of short- and 
long-distance dynamics and a clean distinction 
of model-dependent and model-independent fea- 
tures. 

Yet, a systematic and model- independent anal- 
ysis of exclusive radiative decays is possible in 
the heavy quark limit rrib ^ Aqcd- The rele- 
vant hadronic matrix elements of local operators 
in the weak Hamiltonian simplify in this limit 
because perturbatively calculable hard scattering 
kernels can be separated from nonperturbative 



form factors and universal light-cone distribu- 
tion amplitudes. A power counting in Aqcd / ?7ib 
implies a hierarchy among the possible transi- 
tion mechanisms and allows to identify leading 
and subleading contributions. In particular, ef- 
fects from quark loops are calculable rather than 
being generic, uncalculable long-distance contri- 
butions. Our approach is similar in spirit to 
the treatment of hadronic matrix elements in 
two-body non-leptonic B decays formulated by 
Beneke, Buchalla, Neubert, and Sachrajda 
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2. B ^Vj at NLO in QCD 

The effective Hamiltonian for b 
tions reads 

V ^ p=u,c ^ i=l j=3 ^ 



where Ap ' — VpsVph- The relevant operators 
are the current-current operators Q\ 2, the QCD- 
penguin operators Q3...6, and the electro- and 
chromomagnetic penguin operators Q7,8. The ef- 
fective Hamiltonian for & — > ^7 is obtained from 
(|^) by the replacement s d. The most difficult 
step in computing the B — > decay amplitudes 
is the evaluation of the hadronic matrix elements 
of the operators in (|^). In the heavy-quark limit 
a systematic treatment of the hadronic matrix el- 
ements is possible. In this case the following fac- 
torization formula M is valid 
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Figure 1. 0{as) and leading power contri- 
bution to the hard-scattering kernels T/^ from 
four-quark operators Qi (left) and from Qg. The 
crosses indicate the places where the emitted pho- 
ton can be attached. 



where e is the photon polarization 4-vector. Here 
pB~*v is a, S ^ y transition form factor, and 
^B, are leading twist hght-cone distribution 
amplitudes of the B meson and the vector meson 
V, respectively. These quantities are universal, 
nonperturbative objects. They describe the long- 
distance dynamics of the matrix elements, which 
is factorized from the perturbative, short-distance 
interactions expressed in the hard-scattering ker- 
nels T/ and T^". To leading order in QCD and 
leading power in the heavy-quark limit, Qj gives 
the only contribution to the B —^ amplitude. 
At 0(as) the operators (3i...6 and Qg start con- 
tributing and the factorization formula becomes 
nontrivial. 

The relevant diagrams for the NLO hard- vertex 
corrections T/ were computed in to get the vir- 
tual corrections to the matrix elements for the in- 
clusive 6 — > S7 mode at next-to-leading order. We 
re-interpret these results as the perturbative type 

I hard-scattering kernels for the exclusive process. 
As required for the consistency of the factoriza- 
tion formula, these contributions are dominated 
by hard scales to;, and are hence infrared finite. 

We now turn to the mechanism where the spec- 
tator participates in the hard scattering. The 
non-vanishing contributions to T/^ are shown in 
Fig. |l|. We can express both the type I and type 

II contributions to the matrix elements {Qi) in 
terms of the matrix element (Qt), an explicit fac- 
tor Us, and hard-scattering functions Gi and Hi 
which are given explicitely in |6j. 
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Figure 2. Dependence of the branching fractions 
B{B^ K*°j) and B{B- p'j) on the renor- 
malization scale fi. The dotted line shows the LO, 
the dash-dotted line the NLO result including 
type-I corrections only and the solid line shows 
the complete NLO result. 
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V"f amplitude then can be writ- 
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where the factorization coefficients 07(^7) con- 
sist of the Wilson coefficient C7 and the contribu- 
tions from the type I and type II hard-scattering 
corrections. We get a sizeable enhancement of 
the leading order value, dominated by the T^- 
type correction. The net enhancement of 07 at 
NLO leads to a corresponding enhancement of the 
branching ratios, for fixed value of the form fac- 
tor. This is illustrated in Fig. |^, where we show 
the residual scale dependence for B{B K*'^'-f) 
and B{B~ — > p~7) at leading and next-to-leading 
order. 

A complex phase is generated at NLO which 
leads to a non-vanishing value of the CP asym- 
metry 



AcpiV-/) = 



r{B V-f) - T{B V-f) 
r{B Vj) + T{B Vj) 



It is of 0(10%) for the pj mode with the 
largest theoretical uncertainty coming from the 
scale dependence. The CP asymmetry is with 
Acp{K*-f) « -0.3% very small for the b ^ sj 
transition because of the large CKM hierarchy 
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Figure 3. The isospin-breaking asymmetries 
A(X*7) and A(p7) as a function of the CKM 
angle 7 with and without the inclusion of QCD 
penguin operator effects. 

A further interesting observable is the charge 
averaged isospin breaking ratio 



A{Vj) 



+ ■ 
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where w = 1 for V" = K* and w 1/2 for F = 
p. Within our approximations, isospin breaking 
is generated by weak annihilation contributions. 
They are power suppressed but can neverthe- 
less be computed within QCD factorization be- 
cause the colour-transparency argument applies 
to the emitted, highly energetic vector meson in 
the heavy-quark limit. Isospin breaking was al- 
ready discussed in partially including NLO 
corrections. Kagan and Neubert found a large ef- 
fect from the penguin operator Qq on the isospin 
asymmetry A{K*^) Q. This brings the predic- 



tion A{K*j) = (-7.5+5:9)% within the large er- 
ror bars in rather good agreement with the ex- 
perimental value A{K*"f)exp = (-19.2 ± 11.8)% 
|6|. For B — !■ p7 we find a strong dependence 
of the isospin asymmetry on the angle 7 of the 
unitarity triangle. This is illustrated in Fig. ^ 
A measurement of the isospin asymmetry A{pj) 
can therefore give valuable indirect information 
on the CKM angle 7. 



3. The Radiative Decays B 
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The double radiative B — s- 77 modes realize the 
exceptional situation of nontrivial QCD dynamics 
related to the decaying B in combination with 
a completely nonhadronic final state and simple 
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Figure 4. The sublcading power IPI diagram. 

two-body kinematics. As the two-photon system 
can be in a CP-even or CP-odd state we can study 
direct CP-violating effects. 

In the heavy-quark limit we get a factorization 
formula |^,^ for the hadronic matrix elements of 
the operators in the effective Hamiltonian (|l|): 

(7(ei)7(e2)|Q.|5> = f' d^Tr{0'^>Bi{0^i^e2. 
Jo 

Because there are no hadrons in the final state 
only one type of hard-scattering kernel T (type II 
or hard-spectator contribution) enters the factor- 
ization formula. 

To leading power in Aqc d / "m-b and in the lead- 
ing logarithmic approximation of QCD only one 
diagram contributes to the amplitude for B 
77. It is the one-particle reducible (IPR) dia- 
gram with the electromagnetic penguin operator 
Qt, where the second photon is emitted from the 
s-quark line. The result for both the CP-even and 
CP-odd amplitude is 

A± = -(Ai^)+A(''))C7^ 

where = 0{Aqcd) parametrizes the first neg- 
ative moment of the B meson wave function. In 
the present approximation the strong-interaction 
matrix elements multiplying aI'''' and A^'^ are 
identical and have no relative phase. For the 
strictly leading-power result we therefore get no 
direct CP violation. 

Subleading contributions come for example 
from the IPR diagram, where the second pho- 
ton is emitted from the b quark line. This and 
other subleading effects in (Qy), however, con- 
tribute equally to the up- and charm-quark com- 
ponents of the amplitude and therefore do not 
affect direct CP violation. The one-particle irre- 
ducible diagrams (IPI) of Fig. ^, on the other 
hand, provide the basic effects responsible for 
a difference between the up- and charm-quark 
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sectors of the amplitude, including rescattering 
phases. Although power-suppressed they are cal- 
culable in QCD factorization. We then get a non- 
vanishing CP asymmetry already at ©(a^). The 
effect for Bs 77 is again negligible whereas 
the direct CP asymmetry for Bd — > 77 is ap- 
proximately — 10%. However, we do not expect 
these CP asymmetries to be measured in the 
near future because of the tiny branching ratios. 
We predict B{B, 77) = (1.2^^;^) • lO"*' and 

BiBd ^ 77) = (S.lla.i) ■ 10"* which are both 
roughly two orders of magnitude below the cur- 
rent experimental upper limits. The large uncer- 
tainty in our prediction comes from the poorly 
known hadronic parameter A^. 



4. Conclusions 

We have presented a systematic and model- 
independent framework for the rare radiative de- 
cays B Vj and _B ^ 77 based on the 
heavy-quark limit mi, ^ ^qcd ■ Quark-loop con- 
tributions are calculable in QCD factorization 
rather than being uncalculable long-distance ef- 
fects. Non-factorizable long-distance corrections 
may still exist, but they are power-suppressed. 
Strong interaction phases from both hard-vertex 
and hard-spectator contributions are calculable 
and important for CP-violating observables. We 
have seen that weak-annihilation amplitudes in 
B — > are power-suppressed but numerically 
enhanced and calculable. We used them to esti- 
mate isospin breaking effects. 

Our NLO predictions for the central values 
of the B branching ratios are B{B^ 

K*Ory) ^ {7Atli) ■ 10-^ and B{B- p-7) = 



(1.6_Q 5) • 10~ . They are substantially larger 
than the leading logarithmic values when the 
same form factors are used. The dominant uncer- 
tainty comes from the variation of the nonpertur- 
bative input parameters, most notably from the 
B ^ V form factors. This situation, however, 
can be systematically improved. In particular, 
our approach allows for a consistent perturbative 
matching of the form factor to the short-distance 
part of the amplitude. 
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